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We have calculated the fission probabilities for 237 Np, 233,235,238^ 232 Th, and na t pb following 
the absorption of photons with energies from 68 MeV to 3.77 GeV using the RELDIS Monte-Carlo 
code. This code implements the cascade-evaporation-fission model of intermediate-energy photonu- 
clear reactions. It includes multiparticle production in photoreactions on intranuclear nucleons, 
pre-equilibrium emission, and the statistical decay of excited residual nuclei via competition of 
evaporation, fission, and multifragmentation processes. The calculations show that in the GeV en- 
ergy region the fission process is not solely responsible for the entire total photoabsorption cross 
section, even for the actinides. The fission probabilities are 80 — 95% for 233 U, 235 U, and 237 Np, 
70 - 80% for 238 U, and only 55 — 70% for 232 Th. This is because certain residual nuclei that are 
created by deep photospallation at GeV photon energies have relatively low fission probabilities. 
Using the recent experimental data on photofission cross sections for 237 Np and 233 > 235 . 238 U from 
the Saskatchewan and Jefferson Laboratories and our calculated fission probabilities, we infer the 
total photoabsorption cross sections for these four nuclei. The resulting cross sections per nucleon 
agree in shape and in magnitude with each other. However, disagreement in magnitude with total- 
photoabsorption cross-section data from previous measurements for nuclei from C to Pb calls into 
question the concept of a "Universal Curve" for the photoabsorption cross section per nucleon for 
all nuclei. 

PACS numbers: 25.20.-x, 25.85.Jg, 24.10.Lx 



I. INTRODUCTION 

A. Relevant Framework 

Recently, high-precision experimental data have been 
obtained for the photofission of actinide and preac- 
tinide nuclei at the Saskatchewan Accelerator Laboratory 
(SAL) and Jefferson Laboratory (JLab) HUM- These 
data constitute both a challenge and an opportunity: a 
challenge to theory to reproduce them, and an opportu- 
nity to use them to determine the total photoabsorption 
cross sections for these nuclei. From these one can throw 
light, for example, on the concept of a "Universal Curve" 
for the photoabsorption cross section per nucleon for all 
nuclei pg. 

At this time, the only theory capable of reproduc- 
ing these data is an extension of the widely-known In- 
tranuclear Cascade (INC) model to GeV photon- 
induced reactions 0], including subsequent evaporation 
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and fission. Almost three decades ago, the INC model 
was successful in describing experimental data obtained 
with bremsstrahlung photons Later, the hybrid 

precompound-evaporation model approach of Ref . j9[ was 
successful in describing multiple photoneutron reactions 
below the pion-production threshold. 

A related subject is the electromagnetic dissociation 
of relativistic heavy ions, where virtual photons are used 
to initiate photonuclear reactions [l£J, El • Recently, the 
INC model |7| has been applied to the electromagnetic 
dissociation of ultrarelativistic heavy ions 
In such reactions, nuclei are disintegrated by virtual pho- 
tons over the wide energy range from a few MeV to a few 
tens of GeV. A Monte-Carlo code called RELDIS (Rel- 
ativistic ELectromagnetic Dissociation) was devised to 
perform calculations for real and virtual photons across 
this wide energy range. 

We now use the RELDIS code to describe the re- 
cent photofission data 0, 0, 0]- This paper reports 
our results and their implications. The description of 
the RELDIS model of photonuclear reactions is given in 
Sec. [H] Calculational results for nuclear fission proba- 
bilities and absolute photofission cross sections are pre- 
sented in Sec. IIIII Evaluation of the photoabsorption 
cross section per bound nucleon, commonly referred to 
as the "Universal Curve," from the actinide photofission 
data and calculated fission probabilities is presented in 
Sec. IIIII We summarize our findings in Sec. IIVI 



2 



B. Experimental Background 

As first pointed out by Bohr and Wheeler 0| an d 
echoed by Aage Bohr 0]> the relative simplicity and di- 
rectness of the electromagnetic interaction is very useful 
in the exploration of the process of nuclear fission. At the 
first stage of a photonuclear reaction, the absorption of a 
photon on a pair of intranuclear nucleons brings in only 
thermal energy and produces fewer changes in the struc- 
ture of the target nucleus than, for example, reactions in- 
duced by protons or antiprotons. In the latter cases, the 
antiproton annihilation removes an intranuclear nucleon, 
or, respectively, the projectile nucleon can be trapped 
by the nuclear potential. Therefore, one can expect that 
following photoabsorption at low energies, most of the 
fissioning nuclides are similar to the target nucleus. This 
makes photofission studies much more transparent. 

Over the years, photon-induced fission has attracted 
attention, and considerable progress has been made in 
experimental photofission studies. In the Giant Dipolc 
Resonance (GDR) region, experiments at Livermore [Tsl 
TTgL |20| and elsewhere delineated the giant-resonance pa- 
rameters with high accuracy for eight actinide isotopes 
( 232 Th, 237 Np, 239 Pu, and the five long-lived uranium iso- 
topes). At these low energies (from threshold to about 
20 MeV), photofission reveals its simplicity, and the to- 
tal photoabsorption cross section equals the sum of the 
single- and double-photoneutron cross sections plus the 
photofission cross section, since the high Coulomb barrier 
greatly inhibits the emission of charged particles. 

In the quasideuteron and A33(1232) regions, experi- 
ments at Saclay [2lJ , Mainz , and elsewhere extended 
the photofission data to intermediate energies for 2 35,238tj 
and 232 Th. From a comparison of the photofission cross 
sections per nucleon for these nuclei with the total pho- 
toabsorption cross sections per nucleon obtained by other 
methods at Mainz |2^|, Frascati j^, and Bonn [25| . 
it was commonly believed that there was a "Universal 
Curve" and that the photofission cross section for these 
nuclei saturates the total cross section. The recent pre- 
cise data from SAL 0, agree well with the earlier 
data from Saclay and Mainz across the upper part of the 
quasideuteron region (about 60 MeV to the photopion 
threshold) and the lower part of the A 33 (1232) region 
(up to about 250 MeV). 

The recent data from JLab 0,0 extend our knowledge 
of the photofission cross sections for 237 Np, 233 > 235 > 238 U, 
232 Th, and nat Pb across the A 33 (1232) region and 
through the higher-resonance region (from 0.20 to 3.77 
GeV). In this energy domain, photonuclear reactions 
seem to be more complicated. New advanced theoreti- 
cal models must be employed to account for the greater 
number of reaction channels that are open during the 
first step of the reaction due to meson photoproduction 
on the intranuclear nucleons. Also, the understanding 
of the role of nuclear fission among the many other de- 
cay modes of the excited nucleus requires a well-founded 
theoretical model to describe such decay. 



C. Theoretical Approaches to Describe 
Photofission 

One of the first calculations using the Intranuclear Cas- 
cade (INC) model that includes photon absorption at 
intermediate energies was that of Ref. 0. This calcu- 
lation takes into account the channels of the interac- 
tion with production of only one or two pions, and thus 
can be applied at photon energies only up to 1 GeV. 
It was tested with old experimental data obtained with 
bremsstrahlung photons. Later, the model was success- 
fully applied to early photofission studies (2f| ■ 

INC model predictions || |2?J for the first pre- 
cquilibrium stage of the photofission reaction were also 
used in fission probability calculations for 232 Th and 238 U 
nuclei by introducing the mean-compound-nucleus ap- 
proximation [2{j . According to this method, the ensem- 
ble of excited residual nuclei created after the first step 
of a photonuclear reaction is replaced by a single excited 
nucleus with the average values of neutrons N, protons 
Z, and excitation energy E* . In order to explain the data 
on thorium, one needs to account for its higher nuclear 
transparency relative to uranium [29j. Such a difference 
for nuclear systems of comparable mass was attributed 
to subtle details of the nuclear structure of 232 Th and 
238 U. The nature and origin of such differences still need 
to be explained. The fissility of 238 U was assumed to be 
saturated; that is, its fission probability Wf = Tf/<Jtot 
was assumed to be unity. 

Another phcnomcnological method was proposed in 
Ref. H3| to describe the photofission of 209 Bi, 232 Th, and 
238 U in the quasideuteron region, E 1 = 30 — 140 MeV. 
A special phenomenological factor, which selects the 
quasideuteron absorption that leads to photofission, was 
introduced and was found to be different for 209 Bi than 
for 238 U. 

In Refs. |3lL l32T| the photofissility of actinide nuclei 
at intermediate energies was considered within the mul- 
ticollisional model for photon-induced intranuclear cas- 
cade process and the statistical neutron evaporation and 
fission models for deexcitation of residual nuclei. The fis- 
sion probabilities for both 237 Np and 238 U were found to 
be unsaturated, i.e. very close, but still less than unity. 
However, the authors of Refs. HJIIjI considered their ap- 
proach to be accurate only in a restricted range of photon 
energies, from 0.5 GeV to 1 GeV. This makes difficult di- 
rect comparison with the SAL and JLab data 0, S 0. 3 , 
since the data were obtained over a much wider energy 
range, from 68 MeV to 3.77 GeV. 

Recently, the range of applicability of the INC model 
was extended up to 10 GeV 0. The accuracy of the 
model was improved from the photoncutron thresholds 
up to 2 GeV, and new data obtained with monochro- 
matic photons were used to verify the model predictions. 
The analysis of photofission reactions performed in the 
present paper is based on this approach, which is free of 
the simplifications and assumptions made in Refs. 29, 30] 
and is valid for a much wider region of incoming photon 
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energies than the approach of Refs. [3lll3^|. 
The following questions are addressed here: 

(1) Is it possible to describe the photofission of actinide 
nuclei over such a broad energy range, and if so, 
how good is this description? 

(2) Does the relative simplicity of low-energy photonu- 
clear reactions persist at photon energies above 
1 GeV, where multiple photoproduction of hadrons 
takes place? 

(3) To what extent does the excited residual nucleus, 
created after the first stage of the photonuclear re- 
action, retain the properties of the target nucleus? 

(4) Can the total photoabsorption cross section for 
heavy actinide nuclei be obtained from the total 
photofission cross section? 

(5) Is the concept of the "Universal Curve," which has 
been obtained for A < 208 nuclei, also valid for 
A > 233 nuclei? 

(6) What is the level of calculational and experimental 
uncertainties in the evaluation of the total photoab- 
sorption cross section for heavy actinides? 

In order to answer these questions, we use the model with 
its ma in p arameters extracted from independent stud- 
ies 0, EE H3, 0, H3, HE El of photon-, hadron- 
, and heavy-ion-induced reactions, in contrast to other 
approaches j^E I^E EI E3 i where only photon-induced 
reactions were considered. The present theoretical inves- 
tigation is aimed at estimating the characteristics of the 
first step in the photoabsorption process and determining 
the total fission probability for such a reaction. 



II. THE RELDIS MODEL OF PHOTONUCLEAR 
REACTIONS 

A. Initial Interaction and Intranuclear Cascade of 
Photoproduced Hadrons 

The fast hadrons produced in a primary 7 TV or NN in- 
teraction initiate a cascade of successive hadron-nucleon 
collisions inside the target nucleus during the first, 
nonequilibrium, stage of the photonuclear reaction. The 
duration of this stage T cas can be estimated as the time a 
fast particle needs to cross the nucleus: r cas ~ r Q , where 
r < 10- 22 s. 

The INC model is a numerical method to solve the 
kinetic equation that describes hadron transport in the 
nuclear medium. Calculations of intranuclear cascades 
were performed by using a Monte Carlo technique. The 
nucleus is considered to be a mixture of degenerate Fermi 
gases of neutrons and protons in a spherical potential 
well with a diffuse boundary. By using the effective real 



potentials of nucleons and pions, the influence of intranu- 
clear nucleons on these cascade particles is taken into ac- 
count. The momentum distribution of intranuclear nu- 
cleons is calculated in the local-density approximation 
of the Fermi gas model. The distribution of nuclear 
density is approximated by a set of step-like functions. 
The cross sections of elementary collisions NN — > NN, 
NN -> ttNN, ttN -> nN, ttN -> rnrN, ttNN -> NN,. . . 
in the nuclear medium are assumed to be the same as in 
vacuum except that using the Pauli Principle prohibits 
transition of the cascade nucleons into the states already 
occupied by the intranuclear nucleons. A more detailed 
description of the INC model is given in Ref. [(| . 

By means of the INC model, the process of dissipation 
of the initial photon energy can be investigated in detail. 
Part of this energy is transformed into internal excitation 
of the residual nucleus; the rest is released by fast cascade 
hadrons leaving the nucleus. 

B. Creation of an Excited Compound Nucleus 

At the end of the hadronic cascade, some quasiparticles 
remain in the nuclear Fermi gas as "holes" (iV/i), which 
are knocked-out nucleons, and particles (N p ), which are 
slow cascade nucleons trapped by the nuclear potential. 
In the INC model, the excitation energy of the residual 
nucleus is defined as 

£* = X>?+£^ (!) 

i=l i=l 

where the quasiparticle energies for particles e\ and holes 
t\ are measured with respect to the Fermi energy. The 
numbers of nucleons and protons of the residual nucleus 
are given by the relations 

N c 

A RN = A-J2^ (2) 

i=l 

and 

Z RN = Z-Y j e c i , (3) 
t=i 

where A and Z are the numbers of nucleons and protons, 
respectively, in the target nucleus, and q% and el are the 
baryon number and charge carried away by the cascade 
particle i. 

The moment when all fast particles have left the nu- 
cleus marks the beginning of the period of the establish- 
ment of thermal equilibrium in the residual nucleus with 
a duration of r peq ~ (10 — 100) • r Q . At the end of this 
period, thermal equilibrium is reached and an excited 
compound nucleus has been created. 

Various criteria may be applied to decide whether the 
nuclear system has reached thermal equilibrium. Accord- 
ing to the exciton model |36j | , the system of Fermi parti- 
cles created after the cascade stage is not in equilibrium 
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if the number of quasiparticles N q = N p + Nh is less than 
the equilibrium value N^ q ~ y/2gE*, where g is the den- 
sity of single-particle states. Therefore, in the course of 
the equilibration process, pre-equilibrium particles may 
be emitted. In our calculations, the pre-equilibrium exci- 
ton model |3fi| is used to simulate pre-equilibrium emis- 
sion. 

Relative probabilities of the de-excitation processes are 
determined by the excitation energy E* and by the mass 
Arn and charge Zrn of the residual nucleus formed after 
the establishment of thermal equilibrium. 

Depending on the photon energy E 1 , various processes 
may contribute to the energy deposition. In the follow- 
ing, we consider step-by-step all of these mechanisms, 
with special attention to the amount of energy that is 
transformed into internal excitation of the nucleus. 

Figure ^ shows the average excitation energy E* and 
the fraction of E-y that on average is transformed into 
E* in photoabsorption on nat Pb, 232 Th, and 237 Np. The 
average values of E* per nucleon of the residual nucleus 
(E* /Arn) also are shown in Fig.^ 
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FIG. 1: a) Average excitation energy of the residual nucleus 
(E*) following photoabsorption. b) Ratio of the average exci- 
tation energy of the residual nucleus (E*) to the input photon 
energy E 1 . c) Average values of E* per nucleon of the resid- 
ual nucleus. The RELDIS code results are given by solid, 
dashed, and dotted lines for nat Pb, 232 Th, and 237 Np, respec- 
tively. The value for (E*) deduced from an experiment on the 
photoabsorption for lead [3^| is denoted by the data point. 

As shown by calculations, when a nucleus absorbs a 



photon in the GDR region (6 < E 1 < 30 MeV) its en- 
ergy is almost completely transformed into excitation en- 
ergy E*. For a preactinide nucleus like Pb, whose fis- 
sion threshold is above 20 MeV, the de-excitation pro- 
ceeds mainly through the evaporation of neutrons, since 
their separation energies are only about 7 MeV. Due to 
the high Coulomb barrier in heavy nuclei, proton emis- 
sion is suppressed in the GDR region. Fission thresholds 
are much lower for the actinides; e.g., for U and Np, 
it is < 6 MeV 18|, and the fission channel dominates 
in the decay of such nuclei even at low excitation ener- 
gies mm. 

Starting from E 1 — 30 MeV, where the quasideuteron 
mechanism becomes important, up to the single-pion- 
production thresholds at E~ ~ 140 MeV, only a small 
part of the photon energy is converted (on average) into 
the excitation energy E* of the residual nucleus. The rest 
of the photon energy is taken away by the fast nucleons 
originating from the absorbing pair. It was deduced from 
experimental data [jtfl that E* = 43.4 ± 5 MeV for pho- 
toabsorption on lead at E 1 — 70 MeV. The prediction of 
our model agrees well with this value, as seen in Fig. Qp,. 

The two-nucleon absorption cross section of a photon 
on a heavy nucleus 0^4 is taken from the quasideuteron 
model of Levinger [33 , as modified in Ref . [3jJ : 

a^=kZ(l-Z/A)aT ch . (4) 

Here a^ xch is the meson-exchange part of the cross section 
ad for deuteron photodisintegration, 7c? — > np |4jj, A 
and Z are the mass and charge numbers of the relevant 
nucleus, and k «11 is an empirical constant taken from 
the analysis of Ref. [3^| . 

Although the cross section ad decreases strongly with 
photon energy, the two-nucleon absorption mechanism 
competes with the single-nucleon photoabsorption chan- 
nel 7-ZV — > nN even up to E 1 ~ 500 MeV, when the 
wavelength of the incident photon becomes much smaller 
than the intcrnuclcon spacing. 

An interesting effect concerning the photon energy dis- 
sipation above the jN — > irN threshold was noticed in 
Ref. [ll]. A pion of 50-100 MeV has a small interaction 
cross section with nucleons and therefore has a high prob- 
ability to carry away a large part (~ m w ) of the photon 
energy. Only at E 1 w 250 MeV does the average value 
(E*) start to increase, as shown in Fig. ^ 

The reliability of the model predictions for E* can be 
tested by the comparison of the calculated first and sec- 
ond moments of multiplicity distributions of neutrons, 
(N n ) and W„. = y/ (N%) — (iV„) 2 , with the experimental 
data of Ref. [3jj. In photoabsorption on lead, most of the 
neutrons are emitted via evaporation from excited com- 
pound nuclei, and both (N n ) and Wn are sensitive to E*. 
This test was performed in Ref. 13] and (N n ) and W n 
were found to be described with an accuracy of 10-15% 
at 70 < E 1 < 140 MeV. 

Above the two-pion production threshold, at E 1 ~ 
400 MeV, the photon-nucleon (jN) interaction becomes 
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more complicated because of the opening of many possi- 
ble final states. We use a phenomenological model for the 
description of the 7 V interaction, developed in Ref. Q. 
The model includes not only the excitation of nucleon res- 
onances, but also both the resonance contribution from 
the two-body channels, jN — + ttB* and 7V — > M*N 
(B* and M* being baryon and meson resonances) and 
the nonresonant contribution from the multiple-pion pro- 
duction channels 7 V — > iirN (2 < i < 8). Finally, when 
the photon energy reaches a value of a few GeV, multiple- 
pion production becomes the dominant process. A large 
number (~ 80) of many-body subchannels are included 
in our calculation. 

Since the cross section jN — > hadrons is small com- 
pared to the total NN or ttN cross sections, this leads 
to a large number of photons leaving the nucleus with- 
out producing hadrons in the case of direct use of a 
Monte-Carlo simulation technique. This is a natural con- 
sequence of the well-known fact that nuclei are highly 
transparent to photons (i.e., photons do not interact 
strongly). In order to reduce computation time, we simu- 
late photon absorption in each event, but in this case we 
have to normalize our simulated results to the total pho- 
tonuclear cross section to get the absolute fission cross 
sections. 

At high photon energies, many hadrons are produced 
inside the target nucleus in the course of a cascade pro- 
cess. Due to the knockout of intranuclear nucleons by fast 
nucleons and pions, the slower particles pass through a 
lower-density region, thereby undergoing fewer rescatter- 
ings; this is the so-called "trawling" effect. As was shown 
111 Refs. 0,0 E3> 

the trawling effect is important for 
a realistic description of reactions at projectile energies 
above several GeV. In the present calculation, we use a 
nonlinear version of the INC model , which takes into 
account the local depletion of nuclear density during the 
development of the intranuclear cascade. 

In summary, as can be seen from Fig. ^ when the 
photon energy increases from the GDR region to several 
GeV, the nature of the photoabsorption process evolves 
from the excitation of collective nuclear degrees of free- 
dom to the excitation of a single nucleon inside the nu- 
cleus. In the latter case, up to 95% of the photon energy 
is released in the form of fast particles leaving the nu- 
cleus. Nevertheless, the remaining energy deposited in 
the compound nucleus is sufficient for evaporating many 
neutrons and thus leading to fission, since neutron evap- 
oration increases the fissility parameter Z 2 /A. However, 
Fig. H shows only the average values; the entire distri- 
bution of excitation energy E* is shown in Fig. In 
high-energy photoabsorption the E* distribution is very 
wide, since some of the reaction channels lead to strong 
heating of the nucleus by a multipion system. 

As shown in Fig. the excitation energies of resid- 
ual nuclei produced in photoabsorption on lead, thorium, 
and neptunium arc very close to each other, despite some 
differences in the masses of the target nuclei. Therefore, 
one can expect that the dramatic differences in fission 
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FIG. 2: Distributions of E* and Z%^ n /Arn for residual nuclei 
created after the fast stage of the photoabsorption process 
for E 1 = 68 MeV (top) and for E 1 = 3.77 GeV (bottom), 
for 232 Th, 238 U, and 237 Np (left to right). The same number 
(10 5 ) of residual nuclei is shown for each case. 



probabilities in photoabsorption on such target nuclei are 
due to differences in fissility parameters for these ensem- 
bles of residual nuclei, rather than in their excitation en- 
ergies. Figure |21 confirms this expectation. In Fig. |3 
the values of E* and the liquid-drop-model fissility pa- 
rameter Zftj^/Afiis! are given for the whole ensemble of 
excited residual nuclei created after the cascade stage of 
photoabsorption at E~ values of 68 MeV and 3.77 GeV 
for 232 Th, 238 U, and 237 Np. 

As one can see in Fig. |5J at E 1 = 3.77 GeV, residual 
nuclei are created with very broad distributions in exci- 
tation energy, < E* < 300 MeV, and fissility param- 
eter, 26 < Zft N /AiiN < 38. At higher photon energies, 
these broad distributions are determined mainly by the 
variations of the number of pions and nucleons partici- 
pating in the cascade flow. Such variations are due to 
the wide variety of open channels in the primary 7 V in- 
teraction and in the secondary nN and NN interactions 
with intranuclear nucleons. Taken alone, the variation in 
the number of cascade nucleons due to different impact 
parameters in a photonuclear interaction cannot provide 
such broad distributions. Indeed, the impact-parameter 
distributions of photonuclear interaction events are sim- 
ilar for E 7 = 68 MeV and 3.77 GeV; but in the for- 
mer case, the distributions in Z^jsj/Arn are quite nar- 
row. This is explained by the fact that only a limited 
number of channels are open at E^ = 68 MeV, namely 
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7 + (pn) — > p + n and channels of the NN interaction 
below the pion-photoproduction thresholds. 

At higher photon energy, £L > 1 GeV, the photoab- 
sorption process reveals its complexity and the mean- 
compound-nucleus approximation used in Ref. |2flj | for 
fission-probability calculations breaks down. This is 
demonstrated in Fig. |3| where the fissility-parameter dis- 
tributions of the residual nuclei are shown for E y = 
68 MeV and 3.77 GeV. Strictly speaking, the approxima- 
tion used in Ref. [2^| is not valid even at E 1 = 68 MeV, 
since basically two types of photonuclear reactions oc- 
cur: preferential removal of protons, which decreases 
Zrn I 'Arn ', and preferential removal of neutrons, which 
increases of Z\ N / 'Arm . Therefore, the ensemble of resid- 
ual nuclei created at low energies, particularly below the 
pion-production thresholds, can be represented by a pair 
of "characteristic" compound nuclei. The first group of 
nuclei, whose values of Z% n /Arm are less than that for 
the target nucleus, is created mainly via (7,p), (j,pn), 
and (7,p2n) reactions, while the second group is created 
via (7, n), (7, 2n), and (7, 3n) reactions. However, as 
shown in Fig. |3| this approximation for the ensemble of 
residual nuclei becomes invalid in the GeV region, where 
a very wide set of residual nuclei is created following the 
fast stage of photoabsorption. An example of such a dis- 
tribution is shown in Fig.0]for photoabsorption in 237 Np. 
In the course of the cascade process, the target nucleus 
loses up to ~ 30 units in charge and up to ~ 70 nuclc- 
ons. Highly excited preactinide nuclei, like Pb or Au, can 
be created in photoabsorption reactions at high photon 
energies. Poorly explored nuclei with A and Z located 
between Po and Th are also represented in Fig. 0] In 
our calculations, we take into account the whole ensem- 
ble of excited compound nuclei without replacing it by a 
"characteristic" average compound nucleus. 



C. Decay of the Excited Compound Nucleus 

As soon as statistical equilibrium is reached in the 
residual nucleus, the statistical approach for nucleon and 
light-particle evaporation and nuclear fission is an ap- 
propriate scheme for calculation of the relative probabil- 
ities of different decay modes of the compound nucleus. 
Such statistical decay of the compound nucleus is the 
slow stage of the photonuclear reaction, having a charac- 
teristic time r ev f 3> t . Here we use the term "residual 
nucleus" to mean the nuclear residue formed at the end 
of the intranuclear cascade, while the term "compound 
nucleus" is used to describe a residual nucleus after the 
establishment of thermal equilibrium. 

In the domain of modest excitation energies, which are 
typical for heavy nuclei, with E* /Acn < 2 — 3 MeV, one 
can include the competition between evaporation and fis- 
sion only; the onset of multifragmcntation is located well 
above this range at E* /Acn > 4 MeV [35J. As can be 
seen in Fig. El there are no residual nuclei with E* /Arm 
exceeding ~ 1.3 MeV, even at E 1 = 3.77 GeV. Therefore, 
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FIG. 3: Probability distributions of Z\ n /Arn for residual 
nuclei created after the fast stage of the photoabsorption pro- 
cess, for £ 7 = 68 MeV (left) and for E 1 = 3.77 GeV (right), 
for 237 Np, 238 U, and 232 Th, top to bottom. The values of 
Z 2 1 A for the target nuclei are shown by the arrows. 



we do not expect any contribution of the multifragmen- 
tation process to photoabsorption on heavy nuclei. 

According to the standard Weisskopf evaporation 
scheme the partial width Tj for the evaporation of 
a particle j — n,p, 2 H, 3 H, 3 He, or 4 He is given by 



_ (2 Sj + l)nj 
3 7t 2 pcn(E*) 



aj nv (E) Pj (E* - B j - E)EdE 



where Sj, (j,j, Vj, and Bj are the spin, reduced mass, 
Coulomb barrier, and binding energy of the particle j, 
respectively. cr J inv (E) is the cross section for the inverse 
reaction of the capture of the particle i to create the 
compound nucleus, pcn and Pj are the nuclear level 
densities for the initial and final (after the emission of 
the particle j) nuclei, respectively. Natural units, with 
h = c = 1, are used in this paper. 

The Bohr- Wheeler statistical approach [16| is used to 
calculate the fission width of the excited compound nu- 
cleus. This width is proportional to the nuclear level 
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FIG. 4: Distribution of masses and charges of excited resid- 
ual nuclei created after the fast stage of the photoabsorption 
process for E 1 = 3.77 GeV in 237 Np. 



density pf at the fission saddle point: 



E —Bt 



2npcN 



(6) 



where Bf is the fission barrier height. Masses and ki- 
netic energies of fission products are calculated based on 
corresponding approximations |43j to experimental data 
which describe the transition from the asymmetric fission 
mode to the symmetric one, along with other features 
of the fission process. Since the masses of fission frag- 
ments were not measured in the experiments reported in 
Refs. we do not consider such distributions in 

this work. 

The decay of the excited compound nucleus is simu- 
lated using the Monte-Carlo method. The competition 
between the various decay channels at each step of the 
evaporation chain is determined by the relation between 
the partial widths for particle evaporation and fission, 
Eqs. JSJ) and JfjJ), respectively. Finally, in order to cal- 
culate the fission probability Wf, the total number of 
fission events in a computer run is counted and divided 
by the total number of simulated photoabsorption events. 
Evaporation from excited fission fragments is also taken 
into account and was found to be negligible. 

In the present paper, we also take into account the mi- 
croscopic effects of nuclear structure in the nuclear-mass 
and level-density formulas, according to Refs. 0,0,0]. 
Such effects reveal themselves as a noticeable difference, 
up to ~ 10 — 15 MeV for heavy closed-shell nuclei, be- 
tween the values of the measured nuclear ground-state 



masses M exp (Z,N) and those predicted by the macro- 
scopic liquid-drop model Mld(Z, N) 0. Moreover, this 
difference in mass SW gs = M exp (Z, N) — M LD (Z, N), the 
so-called shell correction, and the level-density parame- 
ter a (used in Eq. Q below) are strongly correlated. For 
closed-shell nuclei, the actual values of the level-density 
parameter are substantially lower than the average values 
of A/ 8 — A/ 10 MeV -1 , and these values depend strongly 
on the excitation energy. Proper accounting for these ef- 
fects, as well as for pairing effects, is important mainly 
at low excitations, for E* ~ 10 MeV. Although these 
shell effects are very pronounced at low excitation ener- 
gies, they disappear for E* > 30 MeV HHH]. Several 
phenomenological approximations of the level-density pa- 
rameter were proposed in order to account for such be- 
havior. Our calculations are based mainly on the re- 
sults of Ref. ^(|, where the data on the level densi- 
ties, decay widths, and lifetimes of excited nuclei with 
2 < E* < 20 MeV have been analyzed in the framework 
of the statistical model. 

We use the Fermi-gas expression for the nuclear level 
density at excitation energy E*: 



p(E*) 



12E* 5 / i a 1 / i 



exp{2\/ZE 7 } 



(7) 



where a = TT 2 gp/(S is the nuclear level-density parame- 
ter, which is proportional to the density of single-particle 
states qf at the Fermi surface. 

Pairing-energy effects are accounted for by the substi- 
tution E* — ► E* — A, where the pairing energy A is given 
by: 



A = x ■ 11/VA MeV 



(8) 



with x = 0, 1, or 2 for odd-odd, odd-even, or even-even 
nuclei, respectively. 

The nuclear level-density parameter is a function of Z, 
N, and E* of the corresponding nucleus [44l l4fll l4flj: 



a(Z,N ) E*) = a(A)i l + 6W gs (Z,N) 



where 



, /(£*- A) 
E* - A 



f(E) = 1 - exp(- 7 £); 



a{A) = aA + (3A 2/3 B S 



(9) 
(10) 
(11) 



is the asymptotic Fermi-gas value at high excitation ener- 
gies, and SW gs (Z, N) is the shell correction in the nuclear 
mass formula. The coefficients a and (3 correspond to the 
volume and surface components, which, along with 7, are 
taken to be phenomenological constants. We use the val- 
ues a = 0.114, = 0.098, and 7 = 0.051 (all in MeV -1 ), 
correspondin g to the first set of systematics given in Ta- 
ble 3 of Ref. |4Q| , which includes the shell corrections of 
Ref. 01 • is the surface area of the nucleus in units 
of the surface for a sphere of equal volume. B s w 1 
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for nearly spherical nuclei with small deformation. Since 
the surface is systematically larger at the saddle point, 
B s > 1 there; in fact B s w 2 1 / 3 w 1.26 for the config- 
uration corresponding to the splitting of the compound 
nucleus into two equal spherical fragments. Such a case 
nearly corresponds to the fission of low-fissility nuclei, 
while for highly fissile nuclei the shape at the saddle point 
is closer to a sphere and B s «l. 

In Ref. |42j, for example, the level-density parameter 
aj at the fission saddle point (in the transition state) 
is calculated using an analogous parameter for the neu- 
tron emission channel a n . In this way, a constant ratio 
r = cif/a n is assumed and used as a fitting parameter 
of the model. As a consequence of this assumption, the 
shell-effect influence on the level density for the neutron- 
emission channel is transferred to the level density at the 
saddle point. However, contrary to Ref. J""J, we expect 
that the shell corrections at the saddle point with large 
deformation have no relation to those at the ground state 
with relatively small equilibrium deformation [""""], [""""J, """""J. 
We conclude that the shell corrections should be much 
smaller at the saddle point and therefore we can use the 
asymptotic value a/ (A) instead of an energy-dependent 
af(Z,N,E*). The value of (if (A) is assumed to be pro- 
portional to the asymptotic one for the neutron-emission 
channel a n (A). 

The ratios df/d n for some nuclei were calculated in 
Ref. |E3 

based on the liquid-drop model """l| and using 
an expression similar to Eq. , but with an additional 
curvature term. As shown in Fig. "5J the ratios obtained 
from this procedure are very close to unity for the ac- 
tinides (~ 1.02 — 1.04) and are higher for the preactinides 
(~ 1.07— 1.12). Since the value for df/d n used in our 
calculations needs to be valid for a much wider range of 
nuclides, we tried to interpolate between the values tab- 
ulated in Ref. |5J| to obtain the values for other nuclei 
as shown, for example, in Fig. [""J Such a procedure deals 
with the average values of df jd n and neglects possible 
rapid changes of the ratio for individual isotopes. How- 
ever, this is the only method to estimate df/d n for the 
region of nuclei between Po and Th, where experimental 
information on the nuclear level densities is lacking. 

Using this interpolation procedure, we find that the 
photofission cross sections for lead are underestimated 
by theory by a large factor. Therefore, to improve the 
agreement with the experimental data Lj. |j|. w e increased 
the difference between the values of Ref. [5QJ and unity 
by a factor of 1.7 (variant A) or 1.5 (variant B). The re- 
sulting curves are shown in Fig. [""] As one can see, such 
modification is not too far outside the ran ge o f d f / d n val- 
ues which were obtained in Refs. [5C1I521I*"""*] . Moreover, 
the curve of variant A independently obtained from the 
analysis of the photofission data on nat Pb agrees 
very well indeed with the values of Ref. [""""} , the most re- 
cent and comprehensive of these studies of nuclear level 
densities, denoted by curve "3" in Fig. [""] 

Following Refs. [4*3 . *45l l*"""| . the fission barriers are 
computed from the macroscopic liquid-drop component 
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FIG. 5: The closed circles are the 5//a„ ratios calculated in 
Ref. [*""j based on the liquid-drop model and the dot-dashed 
curve marked "1" shows the result of our interpolation of these 
discrete ratios. The two sets, A and £>, of 5//fi n values used in 
our calculations are represented by the solid and short-dashed 
and dotted curves, resp ectively. For comparison, calculational 
results from Refs. J52, <5.3j are shown by the long-dashed and 
dotted curves, marked as "2" and "3," respectively. 



Bld{Z, N) and the shell correction at the nuclear ground 
state 5W gs (Z,N) [47|: 



B f (Z,N,E*) = B LD (Z,N)-SW gs (Z,N). (12) 



In this expression, the shell correction a the saddle point 
is neglected. 

The validity of the above-mentioned assumptions (ex- 
cept for our interpolation procedure for df/d n ) was con- 
firmed in previous studies of nuclear fission induced by 
various projectiles. For example, the dependence of Wf 
on the target-nucleus mass for reactions with stopped tt~ 
mesons, photons, protons, and a-particles with kinetic 
energies below 1 GeV was considered in Ref. [2(| within 
the framework of the intranuclear cascade, fission, and 
evaporation models. The calculated values were com- 
pared with the experimental data available at that time. 
The influences on fission probability of shell effects, pre- 
cquilibrium emission from the residual nucleus, and the 
parameters of the liquid-drop model were studied. Later, 
the Wf values were deduced within this model for heavy 
nuclei (Z 2 /A > 30) in reactions with these projectiles 
and with pions in flight "27*) . 
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D. Results for electromagnetically-induced 
heavy-ion fission 

As shown above, when energetic photons are absorbed 
by heavy nuclei, a wide range of excited residual nuclei 
is created, having masses and charges far from the sta- 
bility line. Beams of relativistic radioactive nuclei pro- 
vide new access to fission studies of such nuclei. For 
example, in Ref. jH^], fission of 58 secondary projectiles 

^231-234tj 226-231p a 221-229^ 215-226^, 211-223p^ a 
208-212,217,218^ 205-209^ an d' 205 < 206 At) produced by 

fragmentation of 1 GeV/nucleon 238 U nuclei has been 
studied. The total fission cross section in secondary 
collisions of these nuclei with a lead target at energy 
420 MeV/nucleon was found to consist of two comparable 
parts, due to hadronic and electromagnetic interactions. 

Electromagnetically-induced fission results from ab- 
sorption of virtual photons. According to the RELDIS 
model, which is based on the Weizsacker- Williams 
method of equivalent quanta, such nuclei are typified by 
low excitations, (E*) ~ 12—13 MeV, due to the pho- 
tons in the GDR region with energies below 20 MeV. 
The recent data [54[ make possible a crucial test of the 
model, since E* ~ Bf in such processes. This makes the 
model results very sensitive to the fission barrier shape 
and height, as well as to the level-density parameteriza- 
tion p(E*) used in the calculations. 

The experimental fission probabilities w ED for each 
of the ions were derived by dividing the measured fis- 
sion cross sections of Ref. |54| by the calculated total 
electromagnetic-dissociation cross sections a ED . The cal- 
culations were performed by accounting for absorption 
of either one or two photons in each collision event; see 
Refs. 0, 0, for details. Double-photon absorp- 
tion processes were taken into account by applying the 
harmonic-oscillator ansatz in conjunction with the fold- 
ing model |55| . 

Since experimental data on the total photoabsorption 
cross sections for unstable nuclei are not available, the 
cross sections for nearby stable nuclei were used in cal- 
culations of <j ed . Approximations of the total pho- 
toabsorption cross sections obtained in Ref. |5(| were 
used. Such a substitution generally leads to 2-3% er- 
ror in <j ED according to the estimates based on the GDR 
sum rule EEEl. 



Experimental and calculated values of wf D are shown 
in Fig. [HI Very good agreement is found for the most fis- 
sile nuclei 231 - 234 {J an d 226_231 Pa. The description of the 
data for Th nuclei is poor. The calculated values for Ac, 
Ra, Fr, Rn, and At are much lower than the experimen- 
tal data, but good agreement is obtained for the values 



of 



.ED 



for those isotopes which are closest to the stable 



isotopes of Ac, Ra, and Fr. This is explained by the fact 
that the input calculational parameters were adjusted to 
describe the fission of nuclei close to the stability line. 
The calculations presented in Ref. [54[ for Ra nuclei also 
underestimate the electromagnetic fission cross sections. 
At the present time, the description of the fission of ra- 



dioactive nuclei seems to present a common difficulty to 
all fission models. Further work on our fission model 
should be aimed at a better description of the fission of 
unstable nuclei. 
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FIG. 6: Fission probabilities in electromagnetic dissociation 
of heavy ions at 420 MeV/nucleon on a lead target. Exper- 
imental fission probabilities extracted from the fission cross 
section of Ref. |54l| are shown by the full triangles. Calcula- 
tional results are given by the open circles. 



Using single-humped fission barriers instead of double- 
humped barriers and neglecting the collective effects in 
the nuclear level-density formula does not lead to dis- 
agreement with the data for U and Pa nuclei. Below we 
estimate how the failure to describe the fission of other 
radioactive nuclei may affect the total calculated fission 
probability Wf in phototoabsorption. The 58 radioactive 
nuclei listed are only a small subset of all of the residual 
nuclei created in photoabsorption of GeV photons. As a 
rule, such residual nuclei are much more highly excited. 
However, photoabsorption at lower energies, such as at 
E 1 = 68 MeV, creates nuclei with (E*) of only about 20 
MeV, much closer to nuclei with (E*) ~ 12 - 13 MeV, 
which is the case for electromagnetically induced heavy- 
ion fission. This makes useful the comparison presented 
in Fig. □ 

The charge distributions of residual nuclei created in 
photoabsorption on 237 Np and 238 IJ at E 1 = 68 MeV are 
shown in Fig.0 together with the calculated fission prob- 
abilities for residual elements. Excited Np, U, Pa, and Th 
nuclei are created in photoabsorption at this energy. The 
calculated and measured values of w ED for such nuclei in 
electromagnetically induced fission are shown in the bot- 
tom part of the same plot. These values are obtained as 
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FIG. 7: Top: charge distributions of residual nuclei created 
in photoabsorption on 237 Np and 238 U at E~, = 68 MeV. Bot- 
tom: fission probabilities for such nuclei. Calculated fission 
probabilities in photoabsorption at E-, = 68 MeV are shown 
by histograms. Calculated and measured fission probabilities 
for heavy ions at 420 MeV/nucleon are shown by open circles 
and filled triangles, respectively. 



average values of the points for individual nuclides, and 
they are different from the values of Wf for photofission, 
due to the difference in E* . As shown, fission probabil- 
ity in electromagnetically induced fission of U and Pa is 
described at the 5% level of accuracy. The agreement for 
Th is poor, but this element is less abundant in the en- 
sembles of residual nuclei; the estimate of the total fission 
probability Wf for the whole ensemble is quite reliable. 



The reliability of the fission model can be quantita- 
tively assessed by calculating the total fission probability 
for the ensembles presented in Fig. [7] In the following two 
hypothetical examples the experimental and theoretical 
values for w^ D were used, but in both cases the residual- 
element abundances were taken for photoabsorption at 
E 1 = 68 MeV. As found for the 7+Np case, W s = 0.98 
for the calculated wf D , compared with Wf — 0.97 for 
the value extracted from experiment. For the 7+U case, 
Wf = 0.93 and Wf — 0.91, respectively, for the two sets 
of input Wf D values. Therefore, the failure to describe 
the fissility of Th leads to an error in Wf of only ^3%. 
In this way, the reliability of the fission model is demon- 
strated for the most fissile and most probable residual 
nuclei created in photoabsorption on 237 Np and 238 U at 
low energies. 



III. RESULTS AND DISCUSSION 



Fission Probabilities for 237 Np, 233,235,238^ 
232 Th, and nat Pb 



As a rule, the liquid-drop model predicts proton-rich 
nuclei to have higher probability to undergo fission [4H 
IBlT ] . The probability w / for a residual nucleus with given 
mass Arm and charge Z^n to undergo fission during the 
last stage of the reaction is shown in Fig. [S] for E 1 = 
3.77 GeV for 237 Np, 238 U, and 232 Th. 
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FIG. 8: Fission probabilities wj of each of the residual nuclei 
created after the fast stage of the photoabsorption process for 
E 1 = 3.77 GeV on 237 Np, 238 U, and 232 Th (top, middle, and 
bottom panels, respectively). 

The probability Wf is defined for each of the nuclides 
created after the INC stage of photoabsorption, while 
the fission probability Wf refers to the appropriately 
weighted average value calculated over the whole ensem- 
ble of residual nuclei. For example, Wf turns out to be 
below 0.2 for some regions of j4_rat and Zrn far from the 
A and Z of the initial target nucleus. However, these re- 
gions do not contribute much to the resulting Wf, since 
the probability to create a nucleus with such Arn and 
Zrn is low, as can be seen in Figs. 0] and i Asa re- 
sult, the Wf values, shown in Fig. |^1 as functions of J5 7 
for 237 Np, 238 U, and 232 Th, are generally above 0.6 de- 
spite the fact that some of the residual nuclei have a low 
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probability Wf to undergo fission. As is also shown in 
Fig. El this is true for 233 U and 235 U as well, for which 
0.8 < Wf < 0.95, almost as large as the corresponding 
values for 237 Np. 
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FIG. 9: Fission probabilities Wf for phot oabsorpt ion in 
237 Np, 233 U, 235 U, 238 U, 232 Th, and nat Pb. Our calculated re- 
sults for variants A and B are given by the solid and dashed 
lines, respectively; variant A without pre-equilibrium emis- 
sion is represented by the dotted lines. Earlier calculational 
results for 237 Np, 238 U, and 232 Th jjjj are given by thin dot- 



dashed lines; for 237 Np, 238 U Hjj, and for 232 Th [5y| are given 
by thick dot-dashed lines. The circle represents the calcula- 
tional result of Ref. |60|. Values inferred in Ref. (isj from 
the old experimental data of Refs. pH are given by the 
triangles and stars, respectively. 

Our calculated values for Wf for nat Pb are much lower, 
~ 0.10 — 0.15, as can be seen in Fig.EI The values for Wf 
are first obtained for each of the most abundant lead iso- 
topes, 206 Pb (24.1%), 207 Pb (22.1%), and 208 Pb (52.4%), 
and then are averaged with appropriate weights to obtain 
the fission probability for nat Pb. 

As can be seen in Fig. two extreme cases are repre- 
sented by photoabsorption on the highly fissile actinidcs 
237 Np, 233 U, and 235 U and on the pre-actinide nat Pb. At 
E 1 = 68 MeV, fission clearly dominates for the resid- 
ual nuclei created after photoabsorption on 237 Np, 233 U, 
and 235 U. However, the fission contribution gradually 
decreases with increasing Ej, since, on average, a wider 
distribution of residual nuclei created at higher E 7 has 
lower fission probability and nucleon evaporation be- 
comes more important. A completely different tendency 
is found for the photofission of lead, which has a clear 
threshold behavior. Due to high fission barriers for the 
residual nuclei created after photoabsorption on nat Pb, 



and only evaporation from excited nuclei can take place 
there. Above E 7 ~ 100 MeV, the fission probability Wf 
gradually increases to Wf ~ 0.1, but this still represents 
only a small part of the deexcitation process. An inter- 
mediate tendency is found in photoabsorption on 238 U 
and 232 Th. For these cases, Wf has a very broad maxi- 
mum at E 7 ^ 0.5 — 1 GeV and decreases gradually with 
Ej above ~ 1 GeV. 

Different mechanisms of photoabsorption play a role 
at different photon energies. As a result, a variation in 
the trend of the average excitation energy versus E 7 is 
found for 180-250 MeV, as seen in Fig.^ The variations 
of Wf{Ery) which are seen in Fig.[5]have the same origin. 
The probability for a pion produced in the 7^ — * ttN 
process to be absorbed in the nuclear medium increases 
rapidly with E 1 as one approaches the region of the 
A 33 (1232) resonance, E 1 ~ 200 - 500 MeV. If this pion 
is re-absorbed in the nucleus, the excitation energy in- 
creases and, therefore, so does the fission probability. 

If the pre-equilibrium emission process after the cas- 
cade stage of photonuclear reaction is neglected in the 
calculations, the decay of the excited compound nucleus 
takes place with higher excitation energy E* and fissility 
parameter Z 2 /A. This generally leads to higher values of 
Wj, as shown in Fig. by dotted lines for 237 Np, 238 U, 
232 Th, and nat Pb. 

The changes in Wf due to the variation oiat/a n {A 
and B) are of the order of 5% for 237 Np, 233 U, 235 U, and 
238 U and 10% for 232 Th, as also can be seen in Fig. ED 
This important observation of the relative stability of the 
calculated Wf is true for actinide nuclei only. This makes 
it possible to infer the total photoabsorption cross section 
from the photofission data |lj, y, |^ and the calculated 
values for Wf. 

Fission probabilities estimated in Refs. pll |49| from 
previous experimental data on absolute photofission cross 
sections are given in Fig.[5]for comparison, although they 
are in poor agreement with each other. The calculated 
Wf value from Ref. |60| is in good agreement with our 
results for 238 U. 

Compared with the approach of Refs. 0, E3- sev- 
eral simplifications were adopted in our fission calcula- 
tions. For example, we do not use double-humped fis- 
sion barriers for actinides, nor do we take into account 
any collective effects in the nuclear level-density formula 
because (i) information (either experimental or theoret- 
ical) on such effects is not available for the whole range 
of residual nuclei, and (ii) at high excitation energies 
(E* > 30 MeV), which are our interest here, the fission 
probability Wf is sensitive mainly to the ratio a//a„. 
The uncertainty in the choice of a//a n exceeds the uncer- 
tainties resulting from the selection of other calculational 
parameters describing the decay of the excited residual 
nucleus. 

As shown in Fig. the difference between our ap- 
proach and that of Refs. [4f| |58| does not lead to a dra- 
matic difference in Wf for 237 Np and 238 U. In the limited 



the fission process is suppressed below E y ~ 100 MeV energy region 60 < E~ t < 240 MeV, Wf ~ 85 — 95% and 
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~ 70 — 85%, respectively, for both approaches. However, 
for 250 < E 1 < 1200 MeV, our fission probabilities for 
232 Th (variants A and B) are noticeably lower than those 
of Ref. [dSj, where the effect of pre-equilibrium emission, 
the trawling effect, and some multiplc-pion photoproduc- 
tion channels were neglected. 

Our calculated fission probabilities relative to 237 Np 
are shown in Fig. HUI for the uranium isotopes, thorium, 
and lead as a function of E 1 , compared with the experi- 
mental data from SAL 0,0 and JLab Again, only 
minor changes are found due to the choice of a//a„ in 
variants A and B for actinides, whose energy dependence 
of relative fission probabilities is flat. 



233 



U 



: 235 u 

1 - -ir-"- 



238 



1 ; 
0.5 : 



u 



: 232, 



Th 



0.5 - o O D0 °°°! 




FIG. 10: Fission probabilities relative to Np for photoab- 
sorption in 233 U, 235 U 238 U, 232 Th, and nat Pb. Our calculated 
results for variants A and B are given by the solid and dashed 
lines, respectively. The open and closed circles represent the 
data from SAL 0,0 and JLab 00 respectively. 

A better description of the experimental fission prob- 
abilities relative to 237 Np is found for 233 U, 235 U, and 
238 U. Larger uncertainties in the calculations exist for 
232 Th, and especially for nat Pb, where the Wt values are 
farther from unity. 

As shown in Figs. I§1 and llOl our calculated fission prob- 
abilities for 237 Np and 238 U are in agreement within ~ 
10-20% with the theoretical results of other authors and 
the estimates based on earlier experimental data from 
the literature, which are available mainly at low energies, 
E 1 < 300 MeV. The fission probabilities for 237 Np and 
233,235, 238tj are s t a hle within ~ 7% with respect to the 
choice of calculational parameters. At E~ > 100 MeV, 



the main uncertainty in our results for these nuclei is due 
to the choice of df/a n . 

For the case of 232 Th, we see a larger discrepancy with 
the data than for the other actinides, while the energy de- 
pendence is still reproduced reasonably well. The 232 Th 
nucleus deviates from the pattern common to the heavier 
actinides at lower energies with regard to fission proba- 
bility |l8j and the number of prompt neutrons emitted 
per fission [l^ because it has a more spherical shape 0] . 
It may be that this deviation persists to GeV energies as 
well. 



B. The Total Photoabsorption Cross Section as a 
Calculational Input 

In order to obtain the absolute photofission cross sec- 
tion, the calculated fission probability must be multi- 
plied by the total photoabsorption cross section. In 
the RELDIS code, the values of the total photoabsorp- 
tion cross section are taken from approximations to the 
existing experimental data. In the GDR region, the 
Lorentz-curve fits with parameters from Refs. l56l l57j . 
corrected according to the prescription of Ref. |61| , were 
used for this purpose. Above the GDR region, where 
quasideuteron absorption becomes dominant, the total 
cross section is taken from the estimate of Ref. 39] (based 
on the quasideuteron model |38|.') 

Above the pion-production thresholds, a "universal" 
behavior, o 1 a(E 1 ) cx A, is observed (see Ref. 
for the latest experimental data). This means that the 
total photoabsorption cross section per bound nucleon 
(jyA (Ej ) I A has the same magnitude and the same energy 
dependence for light, medium, and heavy nuclei (C, Al, 
Cu, Sn, and Pb) at least up to E~ t ~ 3 GeV. Therefore, 
having the data for one nucleus, one can calculate the 
cross section for other nuclei. However, in this energy 
region the universal curve cj 1 a[E 1 ) j A is very different 
from the values extrapolated from the cross sections on 
free nucleons, (Za^ ri + Na^„)/A, which are deduced from 
proton Ulill and deuteron data 0]. For E 1 > 3 GeV, 
the universal behavior breaks down, and, for example, 
the ratio ct 1 a{E 1 ) / A for lead is 20-25% lower than for 
carbon due to the nuclear shadowing effect [6(| |6tJ • 
In order to approximate the total photonuclear cross sec- 
tions at E 7 > 3 GeV, we use recent results for <j 1 a{E 1 )/A 
obtained with the Glauber-Gribov approximation within 
the Generalized Vector Dominance Model (see Ref. |65| 
and references therein). 



C. Absolute Photofission Cross Sections for 23T Np, 

233,235,238^ 232^ nat pb 



Our calculated absolute photofission cross sections 

for 237 N P , 233,235,238^ 232^ &nd „atp b are ghown ^ 

Fig. 1 1 1| and compared with the experimental data of 
Refs. US SB. For 237 Np and 233,235,238 U; var i ant s A 
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and B give consistent results. Except for 232 Th, neglect- 
ing pre-equilibrium emission after the INC stage does not 
lead to a large change in the calculated photofission cross 
section; a small increase in the calculated cross section is 
obtained only for E 1 > 1 GeV. When normalized to the 
"Universal Curve" in this way, our calculated absolute 
cross sections agree within ~ 7% with the data for 233 U 
and 235 U, but are underestimated for 237 Np and 238 U 
and overestimated for 232 Th, especially in the A 33 (1232)- 
resonance region, although qualitative agreement for the 
shape of the cross sections is found. It should be stressed, 
however, that these cross sections were calculated using 
"the universal" a 7 A(E-y)/A dependence, which was ob- 
tained for nuclei with A < 208. 



Contrary to the case for the actinides, large uncertain- 
ties in the calculated photofission cross section for nat Pb 
make unreliable any estimate of its total photoabsorption 
cross section from the photofission data @, 0] . But since 
our study is aimed mainly at photofission reactions in 
the actinide region, we did not attempt to obtain better 
agreement with the photofission cross section for nat Pb 
via further adjustment of the calculational parameters 
for each stable isotope of lead or by introducing an em- 
pirical dependence of fission barriers on E* , as was done 
in Refs. [glllilzS- 

D. Evaluation of the "Universal Curve" from the 
Actinide Photofission Data 
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FIG. 11: Absolute photofission cross sections per bound nu- 

Cleon for 237 N P , 233,235,238 U) 232 Th; afld nat pb caku _ 

lated results for variants A and B are given by the solid and 
dashed lines, respectively. Variant A without pre-equilibrium 
emission is represented by the dotted lines. Cross-section 
values given in the plot for nat Pb (only) are shown on a 
separate scale on the right-hand side. The data are from 

Refs. uasa. 



Using the same procedure, we obtain the absolute 
photofission cross section for nat Pb shown in Fig. 1111 on 
a separate scale, since this cross section was found to be 
much smaller than those for the actinides. Although the 
results of the calculation based on set A for the o//a n 
ratios are adjusted to match the data near the peak of 
A33(1232), at higher energies they grossly overestimate 
the photofission cross section. Moreover, neglecting pre- 
equilibrium emission leads to a large increase of the fis- 
sion probability due to the increase of the excitation en- 
ergy of the compound nucleus. This leads to an unrea- 
sonably large enhancement of the fission channel. 



It has been suggested [23, |23|, Ell that one might be 
able to infer the photoabsorption cross section from the 
measured photofission cross section for very heavy nuclei. 
In this manner, the behavior of the "Universal Curve" for 
such nuclei could be investigated. 

We can calculate the total photoabsorption cross sec- 
tion per bound nucleon a-yA(Ey)/A by means of our cal- 
culated values for Wf(E 1 ) and the actinide photofission 
cross-section data of Refs. using 



1 



A 



A 



(13) 



W f (E 7 ) 

The data for each of the target nuclei used in Refs. P, Q , 
except for nat Pb, were interpolated by smooth curves and 
then divided by the calculated Wf. The results of this 
procedure, shown in Fig. 1121 are in agreement within 
- 10% of each other, except for 232 Th. Those for 232 Th 
are obtained with Wf(E 7 ) values that we found to be 
more sensitive to the choice of calculational parameters 
and hence less accurately defined. Reasonable agreement 
is found with the a lA {E~ 1 )/A data of Refs. pi I2I 
at E~ t > 0.7 GeV, but marked disagreement in strength 
(~ 20%) exists in the A33 region. 

In order to understand how the mass and width of 
the A33(1232) resonance are changed in very heavy nu- 
clei compared with other nuclei and the nucleon, we plot 
a 1 A(E 1 )/A for the average of four nuclei, 237 Np, 238 U, 
235 U, and 233 U, as the solid curve in Fig. Altogether with 
the data points for nuclei with A < 208 of Refs. pll^. 
The weighted average total photoabsorption cross section 
for the free nucleon is shown as the dotted curve: 



A 



(14) 



where Z = 92, N = 144, and A — 236 were taken as 
average values for these four nuclei. The total photoab- 
sorption cross sections on the proton er 7P and on the neu- 
tron o\n were taken as polynomial fits of experimental 
data [ij, EH 0| . Comparison of the nuclear data in 
Fig. A] with the nucleon curve demonstrates that the 
A33 is shifted upward in mass slightly and broadened 
considerably. 
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FIG. 12: Total photoabsorption cross section per bound nu- 
cleon for actinide nuclei evaluated from the photofission data 
of Refs. Q,3 and our calculated fission probabilities. Average 
data for C, AL Cu, Sn, and Pb are shown by the open (24[ 
and closed |23| circles. Typical values of the systematic un- 
certainties for the measurements of Ref. [24t are shown for the 
first and last points. Additional experimental data for 238 U 
above 2 GeV |z2 and for Pb below 0.3 GeV |H El El El 
are shown by the open triangles and squares, respectively. 



FIG. 13: Average total photoabsorption cross section per 
bound nucleon for 237 Np, 238 U, 235 IL and 233 U (solid line) 
using the photofission data of Refs. Q,Q and our calculated 
fission probabilities. The dashed line represents the same es- 
timate, but divided by a factor of 1.2. The data points are 
denoted as in Fie . 1121 The total photoabsorption cross section 
on a free nucleon which is obtained as a weighted average of 
a JP EH an d c 7 n HI data is shown as the dotted line. 



IV. CONCLUSIONS 



It is stated in Ref. [24| that the A33(1232)-resonance 
excitation parameters increase with nuclear density due 
to hadronic resonance propagation and interactions; for 
Pb, compared with Li, the A 33 mass Ma increases lin- 
early up to a few tens of MeV, and its width Ta up 
to several tens of MeV. Following this trend, only small 
modifications for Ma and Ta are expected if the mass of 
the target nucleus increases by 30 units, from A = 208 to 
238. When the solid curve is re-scaled to fit the experi- 
mental points close to the maximum, as shown in Fig. 1131 
by the dashed curve, essentially no further modification 
of the mass or the width of the resonance is required. Our 
study shows a far more dramatic effect: the strength of 
the A33(1232) resonance increases by ~ 20%, without 
any noticeable change in its mass or width. 

Finally, one can also see from Figs. ^1 El and 1 131 that 
a shoulder appears for E 1 ~ 0.8-1.0 GeV. This shoulder 
might well reflect the excitation of the iV*(1520) reso- 
nance. Evidence for this iV*(1520) excitation is seen most 
prominently in the structure of the JLab photofission 
cross sections for 233 > 235 U and 237 Np shown in Fig. EJ 



We have calculated, using the RELDIS model, the fis- 
sion probabilities of 237 Np, 233,235,238 U; 232 Th; and „at pb 
as a function of incident photon energy. We have used 
these values, together with the measured photofission 
cross sections from SAL and JLab, to infer the total pho- 
toabsorption cross sections for these nuclei. We now an- 
swer the questions raised in Sec. II CI 

(1) We have been able to describe the photofission 
of the actinide nuclei over a broad energy range. 
The calculated fission probabilities for 237 Np, 238 U, 
235 U, and 233 U were found to be stable within 10% 
with respect to reasonable changes of the input cal- 
culational parameters. 

(2) Our model calculations show that photonuclear re- 
actions become more complex in the GeV energy 
region, when many reaction channels are open. 
Contrary to the naive expectation that for the ac- 
tinides, the higher the photon energy, the closer the 
fission probability Wf should be to unity, our calcu- 
lated fission probabilities for E y > 1 GeV are found 
not to exceed those values for E 1 ~ 50 — 100 MeV; 
they are, in fact, predicted to decrease somewhat. 
At higher photon energies, only a small part of the 
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energy is converted into the internal excitation of 
the absorbing nuclear system. The resulting behav- 
ior of Wf is due to the creation of compound nuclei 
with fissility parameter Z 2 /A much lower than that 
for the target nucleus. 

(3) We show that the residual nuclei after direct and 
pre-equilibrium emission of hadrons resemble the 
target nucleus, in the sense that the fissility pa- 
rameter Z 2 1 A is similar, only for relatively low in- 
cident photon energies {E 1 < 100 MeV). They dif- 
fer markedly for GeV energies. The products of 
photospallation reactions have relatively low prob- 
abilities to undergo fission, but high rates to decay 
via the less collective channels leading to the evap- 
oration of nucleons. 

(4) Our calculations show that for the actinides, fission 
is not the only significant outcome following pho- 
toabsorption at GeV energies. Even at the high- 
est tagged-photon energies (E 7 ~ 4 GeV) used in 
the measurements at JLab, the competition with 
the evaporation process is still important. There- 
fore, we conclude that the photofission cross section 
cannot be substituted for the total photoabsorp- 
tion cross section in the GeV energy region, even 
for those nuclei with the largest fissility parame- 
ters. The photoabsorption cross section must be 
obtained from a calculation of Wf applied to the 
photofission cross section. 

(5) The total photoabsorption cross sections per bound 
nucleon for 237 Np and 233 > 235 > 238 U were inferred us- 
ing the recent photofission data 0, 0] and our cal- 
culated fission probabilities. The resulting "uni- 
versal" curves for (j 7 a(E-{)/A were found to be 
in agreement within ~ 10% of each other. How- 



ever, reasonable agreement with ij 1 a{E 1 ) / A ob- 
tained for nuclei with A < 208 was found only 
for E 1 > 0.8 GeV, while our values for the heavy 
actinides are ~ 20% higher in the region of the 
A33(1232) resonance. 

(6) The calculational uncertainties for fission proba- 
bilities for Np and U are estimated to be at the 
level of 10%, which should be combined with the 
~ 5% accuracy of the measured photofission cross 
sections. Although the disagreement found with 
<j 1 a(E 1 ) I A for A < 208 nuclei cannot readily be 
explained by the presence of these uncertainties, 
one also should take into account the accuracy of 
previous measurements of the "Universal Curve," 
which can be estimated as no better than 5% as 
well. These uncertainties taken together prevent us 
from drawing a definitive conclusion that the "Uni- 
versal Curve," obtained for nuclei with A < 208, 
breaks down for heavy nuclei with A > 233 and for 
E 1 in the A33(1232)-resonance region. However, 
our findings clearly call into question this concept 
and demonstrate the need for direct measurement 
of the total photoabsorption cross sections for the 
heavy actinides. Such measurements would serve 
to confirm or refute the predictions of this paper. 
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